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1 Introduction 

The standard cosmology teaches that the universe is expanding. This implies a steady 
change of the state of matter. Particularly interesting are phase transitions where 
the structure and the relevant degrees of freedom change. Probably the most drastic 
transition was the formation of the early universe. There is a number of ideas how 
did this happen, including the creation from a quantum state (cf. [Q and further 
references there). Since the extrapolation of our quantum field theoretical models to 
high temperature and density points to a maximum temperature of T ~ 10^^ GeV, 
where thermal equilibrium between various particles can be maintained, one can guess 
that the thermal history starts at such a temperature scale and corresponding world 
age of 10^'^'' sec. Immediately afterwards one conjectures that a symmetry breaking 
phase transition happens. At a scale of T ^ 10^^ GeV, or slightly less, the extrapolated 
coupling strengths of the strong, weak and electromagnetic interactions merge, and one 
generally believes that the behavior of matter is to be described within the framework 
of a grand unified theory (GUT). The very structure of GUT is still matter of debate, 
and also the relation to supersymmetric and string theories is not yet settled down. 
Nevertheless, infiation, i.e. a stage of accelerated expansion and a huge blowing off of 
the space is an important scenario, which allows to understand some of the presently 
observed properties of the universe. 

The subsequent evolution appears quite unspectacular for a long time span: expan- 
sion means cooling, and particles with masses in the order of 40% of the corresponding 
temperature disappear due to annihilation. 

The next interesting stage is the symmetry breaking on the electroweak scale, i.e. 
T ~ 100 GeV. Cooling further, at T ~ 160 MeV the quarks and gluons, roaming up to 
this stage freely, become bound in hadrons. This is the confinement transition, which 
might be called the hadrosynthesis stage. 
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Once the temperature falls below 1 MeV the nucleosynthesis starts, creating the 
light primordial elements. In nuclear matter, a liquid-gas phase transition is conjec- 
tured since a long time (cf. and experimental hints have been verified However, 
in the cosmic evolution this transition does not play a role since it is related to large 
baryon density. 

From such a point of view the history of the early universe can be considered as 
sequence of phase transitions. Indeed, based on quite general models, one can develop 
a scheme to combine the dynamics of phase transitions with the cosmic expansion, as 
envisaged in Here we follow, however, a different avenue and consider the specific 
features of individual transitions. 

Among the above mentioned phase transitions the confinement transition from 
quarks and gluons to hadrons has been elaborated by far in most details. This is be- 
cause it is believed to be reproduceable in the laboratory. Indeed, for colliding nuclei 
at sufficiently high energies, the relevant degrees of freedom should be represented by 
partons, i.e. point-like structures. The presently terminating first round of heavy-ion 
experiments at CERN accumulated such a wealth of data which give rise to the inter- 
pretation that in central nuclear collisions a thermalized state of deliberated quarks 
and gluons has been created. Indeed, on February 9, 2000 a press release of CERN 
announced "compelling evidence for the existence of a new state of matter in which 
quarks, instead of being bound up into more complex particles such as protons and 
neutrons, are liberated to roam freely, • • • this state must have existed at about 10 
microseconds after the Big Bang" 

In accordance with this new information wc will focus here on the confinement tran- 
sition after the Big Bang. We will compare schematically the confinement dynamics 
of the Big Bang with that of the Little Bang (i.e. relativistic heavy- ion collisions). By 
using the known details of the transition we demonstrate the interplay of local matter 
properties and global evolution dynamics. More specifically, we analyze how drives 
the local matter state the global expansion, and vice versa, how the expansion causes 
cooling and transformation of matter forms into each other. 

There is some apparent similarity in cosmology and heavy-ion physics, which we 
try to visualize in Fig. ^ From present observations of remote objects in the universe, 
we look somewhat into the past. Combining observational facts, like the distribution 
and the red shift of galaxies, as seen in the left panel of Fig. |l|, one can develop a picture 
of the early stages of the universe. The situation in heavy-ion physics is very similar 
to this. Observing the created hadrons in a very late stage, one tries to extrapolate 
back to the hottest and densest stages. The artist's view of a semi-peripheral collision 
of heavy-ions, displayed in the right panel of Fig. ^ shows a stage where the matter is 
converted into hadrons which still interact for a while before disassembling. However, 
there is a chance in heavy-ion collisions to receive a direct signal of the very early 
stage: the mean free path of real and virtual photons [Q is large in comparison with 
the system size. For such direct probes the problem, however, arises to disentangle 
the primordial signal from a huge background Q. 

The present survey is restricted to a phenomenological framework. That means, 
we do not go into a microscopic theory of phase transitions, rather describing them 
only in terms of the equation of state. We restrict ourselves to the standard model 
and the grand unifying theory models of the strong and electroweak interactions. Only 
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Fig. 1 Expanding matter: The left panel displays a view on a group of remote galaxies as 
seen by the Ifubble telescope (source of the underlying picture: NASA picture gallery in www 
Q), while the right panel shows the distribution of hadrons in a semi-peripheral collision of 
heavy nuclei a few fm/c after the first touch (source: CERN press release in www 



in a few exceptional cases we quote interesting effects occurring beyond the standard 
model of particle physics within supersymmetric or string-theoretical models. 

The paper is organized as follows. In section 2 the phenomenological framework is 
presented resulting in Friedmann's equations for the Big Bang and the corresponding 
equations of matter in the Little Bang. Section 3 deals with quantum chromodynamics 
(QCD), quarks and gluons and their equation of state. Section 4 exercises a simple 
way to construct first-order phase transitions. The basic issues of nucleation theory 
for handling the phase transformation dynamics are presented in section 5. Section 6 
is devoted to details of the cosmic confinement transition. A short remark concerning 
the present status of knowledge on the electroweak transition is made in section 7. In 
section 8 we deal with an attempt of a phenomenological realization of an inflationary 
scenario. The summary can be found in section 9. 

The material represented is based on www scans till end of April 2000. 



2 Phenomenological framework 

Our phenomenological considerations are based on two cornerstones, namely the de- 
scription of 

(i) local matter properties within the framework of thermodynamics and hydrodyna- 
mics, and 

(ii) global cosmic expansion by Einstein's equation and the cosmological principle. 
We are now going to present these issues and define the propositions. 
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2.1 Einstein equations 



The Einstein equations can be written as 




(1) 



K = 



which can be interpreted as gravity, described by geometry or curvature of the space- 
time on the left-hand-side, being caused by matter, which in turn is described by 
the phenomenological energy momentum tensor, Ty, on the right-hand-side. More 
specifically, Rij is the Ricci tensor, R = Rijg^-' the curvature scalar, and gij the 
metric. The Newtonian constant Gn serves as coupling strength of geometry to its 
source, i.e. to matter. In the following we employ mainly units from particle physics, 
i.e. h = c = kB = ^ for the Planck quantum divided by 2tt, and velocity of light, 
and Boltzmann's constant, respectively.^ The energy-momentum tensor can be split 
into a vacuum part, T™'^, and a matter part, T^ , accounting for excitations above the 
vacuum: 



where e™'^ = _pvac gj^gyj-gg ^j^g local Lorentz invariancc of the vacuum. The presently 
observed accelerated expansion of the universe |^ could be caused a dominating vac- 
uum energy e™'^ « 5 x lO^^Og . gm^^ « 0.6 e™* « IQ-^^^M^^ Here, A/pi = 
1.22x10^^ GeV = 2.17x 10"^ g stands for the Planck mass defined by Mpi = ^Tic/Gn. 
The vacuum energy is the missing link to add up all energy forms, including the sub- 
stantial part of dark matter | pT[ , to an amount corresponding to a flat universe with 
e = e'^"''*, which is "predicted" by inflationary scenarios. The present critical energy 
density is /i^l.SS x 10"^^ g cm"^ with h w 0.68 [|l2|. 

From the Friedmann equations (see ( p^ ) below) a less negative "vacuum" pressure 
of p™'^ < — ie™'^ is enough to drive an accelerated expansion. The present vacuum 
energy could be generated by evolving fields, e.g. the tracker fields [|l3|. This would 
cause a noticeable uncertainty in extrapolating back to the early history of the universe, 
because the relation of the vacuum energy density to the energy density of other matter 
forms is rather unsettled. 

If the vacuum energy, sometimes also dubbed quintessence [p^ , would be constant, 
then it can comfortably be included in Einstein's famous cosmological constant A = 
gvac^ — const, and the Einstein equations would read 



In such a case, the vacuum energy can be neglected at early times, i.e. e™'^ ^ e. We do 
not touch upon the questions on the origin of the cosmological constant [l5[ | or why the 
vacuum energy is becoming operative just now, but we assume in the following that 
the present vacuum energy can be neglected at early times in comparison with other 
energy contributions. We mention in advance that in certain cosmic stages appropriate 
vacua are included. 

The conversion constant he = 197.327053 MeV fm is useful for relating length scales and energy 
scales and for translations into other unit systems as well. 



U ~ -'-ij + ^ ij ' 



(2) 



(3) 
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Fig. 2 A contour plot of the temperature distribution in the universe at last rescattering of 
photons as seen by the cosmic background explorer satellite COBE (source: NASA picture 
gallery in www 

2.2 Cosmological principle and Robertson-Walker metric 

The cosmological principle requires homogeneity and isotropy resulting in a simple 
space-time symmetry according to 80(4], or E(3), or S0(3,l) symmetry groups^ and 
the Robertson- Walker metric (e.g. cf. |16| ) 

ds' = dt^ - R{tf (y^^ + '■'^"') ' 

where the only dynamical quantity is the scale factor R{t). The parameter k = 
±1,0 determines whether the universe has a closed, open, or flat geometry. For early 
cosmology and our purposes it is sufficient to consider the case fc = 0. 

The present microwave background radiation corresponds to a blackbody radiation 
with temperature of T = 2.7277 K with relative fluctuations less than 10^'*. The 
cosmic background explorer COBE quantifles these fluctuations, as shown in Fig. ||. 
This COBE map may be understood as a snapshot of the universe at the age of 400,000 
years. While the tiny fluctuations have been amplifled during the last 14 Gigayears to 
the present structures, like galaxies and clusters, the very early stages of the universe 
can be assumed to be very smooth. 



2.3 Perfect fluid 

In the following we describe matter by the energy-momentum tensor 

= eUiUj + QiUj + UiQj + pij, (5) 

where e is the energy density. The stress part, pij, is orthogonal to the four- velocity 
Ui of matter, i.e. pijU^ = 0, as is the energy flow qi. Within this framework, matter 
is described by a fluid thus covering special cases as dust, real fluid, gas, and vacuum. 
In (^) it is assumed that one unique velocity field is common to all subcomponents 

^ We do not consider the Kantowski class which has no local Minkowski limit. 
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of matter. There are two possible gauges of the velocity: either according to Landau, 
with Ui parallel to energy flow, or according to Eckart, with Ui parallel to baryon flow. 
Since we neglect in most examples the baryon charge, the Eckart gauge is applied. 

In the following we mainly consider non-dissipative media, i.e. heat conductivity 
and shear viscosity and possible other gradient terms with non-standard dissipative 



effects (|18[ are not included. The only dissipative effect we consider in subsection 8T 
is the volume viscosity parameterized by the coefficient ^ entering the stress via 

Pi] = {p- S,u',i){-gi.j + UiUj), u'ui = +1. (6) 

The entropy current then reads 

Si = sui (7) 

and has to satisfy s^-i > 0; here s is the entropy density; the semicolon, denotes 
the covariant derivative. (||, ||) with ^ = define the perfect fluid. 

2.4 Equations of motion of perfect fluid in Big Bang 

Discarding any dissipative effect we now present the equations of motion of matter and 
geometry for the Big Bang. Since the cosmic matter consists of several components, 
labeled by a, we take into account 

e = ^ea, _p = ^Pa- (8) 

a a 

2.4.1 Robertson-Walker metric 

The four velocity in comoving coordinates is Ui — 5^ yielding with the metric 

<. = 3|. (9) 

2.4.2 Projection of the contracted Bianclii identity 

The projection of the contracted Bianchi identity UiT^^ -j — yields 

e + 3|(e+p) = (10) 
for the perfect fluid. 

2.4.3 Current conservation 

If there are (3 conserved currents, (n/3U^).i = 0, the conservation of the corresponding 
charges reads 

R 

np + 3—ni3 — npR^ = const. (11) 

R 

The quantities np can represent the baryon density, or the electric charge density etc. 
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2.4.4 Einstein's equations 

For A = the Einstein equations (|l]) read 

R^J - ^.9y-R = 30^4 , C = Mpji (12) 
A consideration of Tn — RToo and the corresponding left hand side of ( p^ yields 



R = y/C'^R'^e-k. (13) 

For fc = these equations can be combined to result in 

R = CRy/e, (14) 
e = -3C(e+p)V^. (15) 

(Notice that for early cosmology the distinction of A; = ±1 or is usually not impor- 
tant.) To integrate the latter equation, an equation of state, p{e), is needed. It should 
be emphasized that the derivation above does not invoke entropy conservation. 
The i = 0, j — component of the Einstein equations ( |l^ ) reads 

R = -^C^R{e + 3p). (16) 

This shows that e + 3p > results in i? < 0, i.e., a decelerated expansion, while 
e + 3p < means i? > 0, i.e., an accelerated expansion. 

2.5 Equations of motion for the Little Bang 

Various dynamical models suggest that in a very high-energetic collision of nuclei 
the initial conditions on proper time hypersurfaces (see left panel of Fig. ^ t = 
y/t'^ — z'^ — const are constant. In suitable coordinates, namely the proper time r and 
the rapidity Y = atanh(z/i), this looks like a homogeneously longitudinally expanding 
fire cylinder (see right panel of Fig. ||) . We neglect the onset of transversal expansion, 
i.e., Vx = Vy = 0. Then the motion is called scale invariant or boost invariant expansion 
[pO[ , which is described by the four-velocity — 7(1,0,0, t;z) with Vz — z/t and 

^=(l_„2)-l/2. 

The corresponding equations of motion for the perfect fluid then read 

e = --(e-|-p) (17) 

T 

from UiT'^ -j = (a dot means here derivative with respect to r), and 

n^r = const (18) 

from {npu^).^i = 0. To integrate ( |l7| ) again the equation of state p{e) is needed. 

Employing a simple equation of state like p = aT** with a ^ 0(1) and the thermo- 
dynamic relations from the next subsection, the dynamical time scales in Little Bang 
and Big Bang can be estimated as 

Little Bang: £i£ ^ r - n - 1 fm/c (19) 
Big Bang: 2+£ ^ ^ iqW f^/^^ (20) 
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Fig. 3 Symmetry of a high-energy heavy-ion coUision. Left paneh Incoming nuclei move 
nearly on the light cone (dashed lines); the thermalized matter has constant energy density 
and baryon density on the hyperboloid of constant proper time r. Right panel: Before the 
onset of transverse expansion the matter expands homogeneously in longitudinal direction. 

where 1 fm/c = 3 x 10~^'^ sec, and we used a temperature scale oi T ^ 0.2 GeV ^ 
2 X 10^^ K. Eq. ( pO[ ) shows that (i) higher energy densities cause faster expansion, and 
(ii) the time scale of the Big Bang is determined by the Planck mass. 



2.6 Thermodynamics 

To complete the list of propositions let us recall the thermodynamic relations for a 
single-component and single-phase medium: 

thermodynamic potential: p{T,fi), (21) 

Gibbs relation: e + p — Ts = /in, (22) 
Euler relations: s = ^, « = (23) 

where /i is the chemical potential and n the related conserved charge density. 

For the early universe, one estimates from present observations and from Big Bang 
nucleosynthesis calculations, ^\t, ~ 10"^" (here ub is the baryon density and 
the photon density), via entropy conservation a baryo-chemical potential ^ ~ 10~^. 
Therefore, for many purposes one can neglect effects of a finite baryon charge, unless 
one is dealing with the evolution of the baryon density itself. The estimates for the 
lepto-chemical potentials are not on such firm grounds due to the non-observability 
of the neutrino background; but one usually assumes similarly small lepto-chemical 
potentials. 

From the quantum statistics of ideal gases one finds the pressure as a thermody- 
namic potential 







/ dx x^ In 


1 ± exp 1 


/o 






(24) 
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Little Bang: 

T 4 




Fig. 4 Left three panels: An artist's view on a central heavy-ion collision (right: the approach 
stage before first touch, middle: highly compressed strongly interacting matter is formed; 
in the center of the blob quarks and gluons are the relevant degrees of freedom, left: the 
system has completely transformed to hadrons and expands. Pictures by S. Bass exploiting 
a transport model). Right panel: A sketch of the temperature evolution till confinement. 



where s labels the particle species, dg their effective degeneracies (being the number 
of degrees of freedom for each boson, and | x the number of degrees of freedom for 
each fermion, respectively); the upper (lower) sign in (|2j) is for fermions (bosons). 
The pressure is dominated by those particle species with masses nig lower than the 
temperature T. This toy model equation of state gives the time evolution of the 
temperature and entropy density 

Little Bang: T = ri(^)^^'\ st ^ const, (26) 
Big Bang: T = (2V3aC^)"^^^ sR^ = const, (27) 

where Ti is the maximum temperature at which the system is thermalized and ti the 
corresponding time. Clearly, the comoving entropies are constant. 

2. 7 Little Bang vs. Big Bang 

One can visualize the similarities and differences of the evolution of matter in the 
Little Bang and the Big Bang as follows. Let us consider in both cases high enough 
temperatures, so that the strongly interacting matter component is in a deconfined 
state. Initially, in the Little Bang (see Fig. ^ the nuclear matter is in its ground state. 
The uni-directional motion during the approach stage has zero entropy. During the 
virulent collision this motion is randomized and a substantial part of the kinetic beam 
energy is converted into particle excitations. According to dynamical models | p9[ , 
after a time scale of ti ~ 0(1 fm/c) or even shorter, the matter is locally equilibrated 
and may be characterized by an initial or maximum temperature Ti . 

According to present estimates the maximum temperatures are 0(200) MeV for 
CERN-SPS| 0(550) MeV for BNL-RHICg and 0(1000) MeV for CERN-LHC| 

Super Proton Synchrotron at CERN, i.e. the European Laboratory for Particle Physics in Geneva 
Relativistic Heavy-Ion Collider at Brookhaven National Laboratory 
Large Hadron Collider at CERN 
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Big Bang: 




Fig. 5 Left panel: A sketch of the expanding universe. Right paneh A sketch of the 
temperature evolution till confinement. 



[p2| . As mentioned above, already at CERN-SPS energies the body of observational 
material is in agreement with the expectation of having produced a state of strongly 
interacting matter where gluons and quarks are the relevant degrees of freedom, ab- 
breviated in the very right panel in Fig. ^ by g, u, d, s. Due to the enormous pressure 
the system expands and converts into hadrons, eventually registered in detectors. 
Note that, despite of the many measured particle species, including the direct probes 
[0, ^, |2^, a simple "yes" or "no" signal, whether a deconfined state is transiently 
created in heavy-ion collisions, is not at disposal. Instead, rather a combination of 
various, partially subtle, observables allows the conclusion to have produced a quark- 
gluon plasma. In contrast to the situation in the cosmic evolution, under laboratory 
conditions the beam energy can be changed, and within some limits the system size 
can be varied as also the impact parameter. 

Presumably, the evolution of matter in the Big Bang (see Fig. |5|) started at much 
higher temperature. Note again the difference of the time scales in Figs. ^ and |^. In 
addition, in the Big Bang there is the non-negligible contribution of the photon and 
lepton background, as indicated in Fig. |^, which is in kinetic and chemical equilibrium 
with the strongly interacting matter. 

There is a simple argument that dense strongly interacting matter must change its 
state. Consider a hadron system at temperature of T ~ 100 MeV. There are mostly 
pions excited with an admixture of kaons, etas, omega and rho mesons, and a few 
nucleons and their anti-particles. At higher temperature, say around T ~ 150 MeV 
the hadrons start touching each other since they have a finite spatial extension on 
the scale of 1 fm (see Fig. ||). At slightly higher temperatures, the hadrons cannot 
be longer the proper degrees of freedom. Instead the matter consists of quarks and 
gluons. Probably also the vacuum structure changes, as indicated by a different gray 
scale in the background of the right panel of Fig. ^. 

Now, we are now interested in the dynamics of the confinement transition, i.e., the 
conversion of deconfined quarks and gluons into hadrons. 
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T = 100 MeV T = 150 MeV T = 200 MeV 

Fig. 6 A sketch of strongly interacting matter at various temperatures. Left and middle 
panels: Hadrons with finite extension are the relevant degrees of freedom. Right panel: Point- 
like quarks and gluons are the relevant degrees of freedom, the vacuum state is changed. 



3 QCD, quarks &: gluons, confinement 

In order to describe the transition from quark-giuon matter to hadron matter in more 
detail one has to rely on the proper theory of strong interaction, i.e. QCD. It rests on 
the Lagrange density (cf. |^ for details, e.g.) 

C^^Piij.D^ ~M)^-\f.^^F^ (28) 

with exact local SU(3) gauge symmetry. is the gauge covariant derivative, V' and 
•ijj stand for Dirac spinors and their adjoints, 7j are Dirac's matrices, and denotes 
the field strength tensor of the non-Abelian gauge fields. Since the different quark 
flavors have different masses, M in ( p8| ) denotes the mass matrix. 

Renormalization on the one-loop level results in the running strong coupling strength 

as(Q') = p7^, f3o = n~ ^Nf (29) 



via dimensional transmutation with the QCD scale Aqcd ~ 0(200) MeV; Nf is the 
number of involved quark flavors. In a system with a very large momentum scale 
Q^, the asymptotic freedom follows immediately from (p9|), i.e., limQ2_,oo — > 0. 
Indeed, if in a gas of quarks and gluons the momentum scale increases with increasing 
temperature, one expects asymptotically an ideal gas. In contrast, at T — s- Tc, the 
system is strongly coupled and perturbative calculations basically fail. 

The only reliable way of ab initio calculations of thermodynamical properties of 
a quark-gluon plasma is to perform Monte Carlo simulations of finite temperature 
QCD in a discretized space-time. A few results of such lattice QCD calculations are 
displayed in Fig. Q together with results of a quasi-particle model |2^ with parameters 
adjusted to the data. 

The transition temperature is estimated as Tc ^ 160 McV ~ ©(m^) ~ 1.6 x 10^^ K 
for a quark-gluon plasma, while for a pure gluon gas it would be Ts'"" ^ 240 MeV (m^ 
is the pion mass). Note that the present lattice QCD calculations cover only the region 
T > Tc and fi = 0. In order to describe the equation of state for hadron matter below 
Tc reliably, a prerequisite would be an accurate description of the mass spectrum of 
the hadrons, a goal not yet fully accomplished. Despite the lack of an equation of state 
at T < Tc, some information on the order of the transition is available. Inspection 
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TIT, TIT, 



Fig. 7 Thermodynamics of a gluon gas (left panel) and a two-flavor (middle panel) and a 
four-flavor (right panel) quark-gluon plasma as a function of the scaled temperature. Lattice 
QCD results (symbols) from the Bielefeld group. The curves represent an adjusted quasi- 
particle model, for details consult [M. 



of the left panel of Fig. || reveals that if the masses of the light current quarks are 
sufficiently small the deconfinement transition is of first order. The physical relevant 
case of the two light fiavors u, d and a medium-heavy s quark is near to the borderline 
to a cross over. The latter notion means a dramatic, but nevertheless smooth change 
of thermodynamic properties in a narrow region around Tc- At sufficiently small 
u, d quark masses and large s quark mass, the phase transition would be of second 
order; the regions of first order and second order are separated by a tri-critical point. 
The pure gauge sector (i.e., infinitely heavy quarks so that they do not longer act as 
dynamical objects) is known to display a first-order phase transition again. 

Since in the early universe the baryo-chemical potential is small the information 
on the region > is not important. Otherwise, in the Little Bang (or more impor- 
tantly, in massive neutron stars with quark cores) the finite baryo-chemical potential 
/J, becomes important. As seen in the right panel of Fig. ||, a richer structure of the 
phase diagram is showing up with di-quark condensates (2SC) and color-flavor locking 
effects (CFL); also an endpoint (E) of the phase borderline may appear at finite values 
of T and ^ (for more details, see [26j). The latter peculiarities make the previous naive 
phase diagram, as shown in Fig. 97mor6 involved. 

To quote some physical numbers characterizing the nuclear matter near ground 
state one can refer to no — 0.15 fm~'^ = 1.5 x lO'^^cm"^, eg = 2.7 x lO^'^gcm"'^, 1 
MeV w IQi" K. 

The deconfinement transition is intimately related to chiral symmetry restoration 
at temperature T^. In hadronic matter, at T < T^, the effective degrees of freedom 
are massive quarks and we have finite expectation values of condensates, like the 
chiral condensate (V'V') ■ For T > the quarks become massless and the condensates 
disappear. For QCD it turns out that deconfinement and chiral symmetry restoration 
coincide, i.e. Tc — T^. Therefore, the chiral condensate may serve as order parameter 
describing the confinement. 
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Fig. 8 Phase diagrams of strongly interacting matter. Left pan el: The order of the phase 
transition as a function of the quark mass parameters (from |25|). Right panel: The phase 
diagram in the T — fj, plane according to fed]. 



Quarks & Gluons 




10 "^"o 



Fig. 9 Naive phase diagram of strongly interacting matter in the T ~ n plane (n is the 
baryon density). 
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4 A constructed phase transition of first order 

In view of the above discussion of the order of the deconfinement transition, it seems 
legitimate to approximate the thermodynamics in the transition region by a suitable 
construction. This construction results in a strong first-order phase transition with 
large latent heat (cf. Fig |ll|) which probably overestimates of what we can expect in 
reality. The construction, however, is made quite general so that it might be applied 
also in other cosmic phase transitions, such as on electroweak and GUT scales to be 
discussed below. A MIT-type bag model is quite popular: asymptotically free quarks 
are confined in a certain spatial region, called bag. The vacuum energy density in the 
bag is enlarged by a constant amount of B when compared with ordinary vacuum. 
The equation of state reads then 

TT^ 4 

Pqcd = dg^u,d-^T — B, dg^u,d = 37, (30) 

with vacuum pressure —B. The use of the pure radiation term cx is suggested 
by asymptotic freedom; near to Tc, the non-perturbative effects are thought to be 
parameterized by B. This ansatz allows only a poor reproduction of the lattice QCD 
results, and for more detailed considerations it must be replaced by other models like 
that in |2^. In some sense this schematic equation of state can be considered as a 
down-extrapolation from asymptotic freedom. Similarly, one can up-extrapolate the 
hadron equation of state, and the most simple ansatz is that of a pion gas neglecting 
the strong interaction: 

P^ = dTr—T'^, = 3. (31) 

yu 

Also this too simple model needs improvements by including the other hadrons (cf. 
[p7|) and their strong interactions. In |Q the interested reader can find a selection 
of some recent approaches and various methods. Despite of many efforts, a generally 
accepted nuclear or hadronic equation of state at large temperature and density has 
not yet emerged. A direct comparison of the various model equations of state with 
experimental data from heavy-ion collisions at BEVALAC, SIS and AGS energies 
(i.e., beam energies of 1---2 and 10 • • • 15 A-GeV) is hampered by dynamical and 
non-equilibrium effects p9| . A recent survey can be found in pO| . 

The equations of state (|3^, ^) are displayed in Fig. |o|. At the crossing of both 
curves, the Gibbs conditions of phase equilibrium are fulfilled, i.e., 

(32) 
(33) 

^^\" = M^^^ (34) 
= Aif, (35) 



Here we consider the special case of /i*^") =0. In the case of competing phases, the 
stable one is that with lower free energy, which is the negative of pressure. The dashed 
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T 



e + 3p<0 

Fig. 10 A constructed phase transition of first order by extrapolating asymptotic equations 
of states. The dashed branches are metastable. 



sections in Fig. are, therefore, metastable branches. According to the Gibbs criteria, 
the conditions for phase equilibrium are fulfilled at Tc — -^jB^^^- 

Fig. |ll| sketches the other thermodynamical state variables as a function of the 
temperature and the energy density. 

Having the two asymptotic branches (|30| , pT[ ) at hand, one can find, of course, 
also smooth interpolations between them displaying a crossover. Or, as exercised 
in ||3^ , one can construct another interpolation with a loop structure in the potential 
p(T); this yields two metastable and one unstable branches, as known from the famous 
Maxwell or double-tangent constructions. Thereby the strong phase transition of first 
order constructed above becomes weaker. However, the mentioned procedures are ad 
hoc and need to be checked carefully against the results of lattice QCD calculations. 



5 Phase transformation dynamics 

In spite of the simplified nature of the phase transition constructed in the previous 
section, we study its implications for the cosmic evolution. As stressed above, this 
construction overestimates the strength of the transition. Therefore, this model serves 
as an upper bound of possible interesting effects caused by the confinement transition. 

Our discussion is based on the classical nucleation theory, where the central quan- 
tity is the critical bubble of the newly forming phase. The change of the Helmholtz 
free energy can be expanded into a power series with respect to the bubble radius r as 

AF{r, T) = — - Ap + 47rcrr^ - 87r7r (36) 
o 

with Ap = Pi—p2 as pressure difference of both phases, a is the surface tension, and 7 
the curvature parameter. Neglecting the latter one, the resulting critical radius, related 
to the maximum of the free Helmholtz energy, is Rn — ; smaller bubbles are 

energetically unfavorable since the surface energy is too large in comparison with the 
volume energy, while larger bubbles are too rare to influence the phase transformation 
dynamics. Critical bubbles are the ones which, after being created by fluctuations, 
can grow further. Notice that the surface tension a is the decisive quantity governing 
the transition dynamics. 
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latent heat 



Fig. 11 Schematic sketch of the thermodynamic state variables for a strong first-order 
phase transition. Upper row: Energy density e and entropy density s as a function of the 
temperature. Lower row: Pressure p and sound velocity = dp/de as a function of the 
energy density. 
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The probability for forming a critical bubble per time and volume unit is given by 

^ = ""'"^{^^^^^1' ^^^^ 
where, from dimensional reasons, 

cr = (TqT^, wq = wqT^ (38) 

must hold, since the relevant scale is given by Tc. One usually assumes that ao,wo ^ 
C(l). The calculation of the coefficient wq is often debated (cf. p3[) , but due to the 
enormously sensitive temperature dependence of the exponent in (|37D, changes in iDq 
may be absorbed in a redefinition of ctq to be discussed below. 

Once a critical bubble is created, it grows with velocity v = v{T, ■ ■ •), where the 
dots indicate a possible dependence on other state variables or transport coefficients. 

Now, we derive an equation for the occupancy of the space with the new phase. 
Thereby, we must take into account that bubbles created at an earlier instant have 
grown and that later a reduced volume for bubble nucleation is available. The general 
growth law x = x(t, w, v, Rc, • ■ •) has therefore also memory effects. Several approxi- 
mations are discussed in pa . Here we take the so-called Avrami approximation 



X2 



l-exp{-/i}, h=— dTw{T)R{Tf 
-J Jo 



r{t) Jr R{e) 



(39) 



where X2 is the volume weight of phase 2. 

Fig. O displays the outcome of a numerical study of the temperature evolution 



during confinement. The equations of state (30, pil) are supplemented by the back- 
ground contribution pbg = "^fagfo-T , with dbg = 14.25. One observes for vanishing 
surface tension a quasi-equilibrium transition with negligible supercooling, see the flat 
section in Fig. ^ Increasing the surface tension parameter ctq causes a deeper su- 
percooling. Then massive nucleation of many critical bubbles happens. This causes 
a sudden reheating to Tc, due to the released latent heat. Then nucleation ceases 
and the transition proceeds via bubble growth. Further increasing the surface tension 
smoothes out the transition and leads to a longer duration of the transition era. In 
the extreme case, hypercooling happens, where the system stays below T^. 

To summarize the typical scales of the confinement transition in the Big Bang let 
us mention that the estimates of the transition temperature Tc ~ 160 MeV from lattice 
QCD yield: 

• beginning of the transition at world age ii ~ 6 /isec, 

• end of the transition in case of a near to equilibrium transition with small surface 
tension t2 ~ 12 /isec. 

Other characteristic quantities are 

• horizon radius Rh ^10 km, 

• Hubble time tn ^ 10^^ sec, 
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Fig. 12 The temperature evolution during a first-order phase transition for various values of 
the surface tension parameter ctq as a function of the scaled dimensionless time r = 2CB^^^t. 
Tl and T2 denote the beginning and end of an equilibrium transition; T° is the temperature 
of the maximum nucleation rate 



• energy density within the horizon Mncorrespondingto ^ IMq, 

• baryon charge within horizon ^ 10'^", 

• scDM ^ 10^^e, ad (from a back-extrapolation of estimates of the present cold 
dark matter energy density). 



6 Details of the cosmic confinement transition 



6.1 Small supercooling 

Early QCD lattice calculations predicted an astonishingly small surface tension ( |3^ ) 
between hadronic and deconfined matter: ctq = 0.0292 |34|, 0.0155 |||, 0.014 • • • 0.03 
(or 0.44) jS^. As these numbers are valid for pure SU(3) gauge theory, where the 
deconfinement transition is undoubtedly of first order, one can speculate whether this 
feature of (Tq ^ 1 is pertinent when the quark degrees of freedom are included. As 
a consequence, the transition is expected to proceed with small initial supercooling, 
AT/Tc < 10^"^, and follows then the unlabeled flat section in Fig. |l^. The temperature 
is kept constant at Tc due to the condensation heat released. Therefore, the transition 
resembles an equilibrium transition since bubbles have time enough to grow and the 
released latent heat is rapidly distributed uniformly. Entropy production is negligible. 
The expansion of the scale factor is slightly changed and points to a mini-inflationary 
era |37|] , i.e., an expansion faster than R oc y/i. 

An opposite scenario is claimed in |38[ . There, within an effective two-phase model 
with strong flrst-order phase transition, the possibility of a substantial supercooling 
is considered and, as a consequence, an exponential growth of the scale factor and 
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Fig. 13 Space-time diagram of two types of growing bubbles of hadronic matter (H) in 
deconfined matter (Q). Left panel: Deflagration bubbles. Right panel: Detonation bubbles. 
The solid lines are the phase boundaries, while the dashed lines depict shock waves reheating 
the matter. The dotted lines indicate stream lines of matter. These sketches are strongly 
stretched in a;-direction so that the light cone is below the dashed lines in (a) and below the 
solid lines in (b). 



a huge entropy increase are advocated. The entropy increase dilutes the baryon-to- 
photon ratio. For instance, a pre-fabricated ratio C'(l) prior to the transition 

can be diluted to the needed value of 10~^ (prior to nucleosynthesis) if the entropy 
increase is tuned to 10^. Scenarios of such a type have been elaborated for inflation 
and will be discussed in section ||. We consider Fig. || and notice the small values of the 
surface tension from lattice QCD. Then the scenario in js^] seems to be too extreme. 
In particular, during the rapid expansion of matter in the Little Bang such parameters 
would cause an extremely deep supercooling, in contrast to the expectation p^ . 

To conclude this subsection, we mention a series of papers [^^, in which the confine- 
ment transition is studied within an effective theory, namely the dual Ginzburg-Landau 
model based on the dual Higgs mechanism by QCD monopole condensation. 



6.2 Turbulent confinement transition? 

There are two types of hydrodynamical solutions for growin g bu bbles in matter being 
at rest far enough from the bubble as depicted in Fig. Deflagration bubbles 
and detonation bubbles. In both types, the bubble growth^ is accompanied by shock 
waves. If we have deflagration bubbles, these waves propagate in the deconfined phase. 
Otherwise, the waves occur in the hadronic phase. As suggested by Fig. ^ there are 
many intersections of the shock waves during the completion of the confinement tran- 
sition. The intersections are related to density inhomogeneities which we will discuss 
below. The density inhomogeneities can represent nuclei for black hole formation visi- 
ble as Hawking flashes or via gravitational lensing or they can cause clumping of dark 
matter. However, the shock waves are very weak and speculations on resulting strong 
turbulences are on less firm grounds. 



See also [M for the bubble dynamics at the end of the cosmological confinement transition. 
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Fig. 14 Schematic view on the process of filling the space with the confined phase by bubble 
growth via deflagration. Many intersections of shock waves arise. 



6.3 Baryon concentration 

Let us now discuss a few possible relics of the cosmic confinement transition. On small 
scales in the order of growing hadron bubbles or vanishing last islands of the quark- 
gluon plasma (cL Fig. |T5|) a baryon concentration arises. According to the scenario of 
a near-to-equilibrium transition, at the phase boundary also chemical equilibrium is 
maintained between quarks (q) and nucleons (N). The reaction <-> iV is unhindered 
and, due to detailed balance, causes the balance equation for the baryo-chemical po- 
tentials 3/ig = UN- Since in the deconfined phase the baryon charge is carried mainly 
by the light current quarks, while in the confined phase the baryon charge is in the 
very heavy nucleons (compared to the temperature scale), the relation tt-qqp 3> 
emerges from the general expression of the net baryon density (density of baryons 
minus density of anti-baryons) in ideal gas approximation (cf. (|2^ , |2^)) 



"■QGP' 



(40) 



where Nf is the number of excited quark flavors (2 or 3 for u, d or u,d,s quarks with 
negligible masses) and Jat = 4 is the nucleon degeneracy; tojv = 939 MeV denotes the 
averaged nucleon mass. ( |40| ) is for /ijv < rriN and in Boltzmann approximation. Due 
to the exponential suppression factor in (^0|), the ratio JT-qcp/^H at T ^ 200 MeV 
becomes very large. In fact, the baryon charge is concentrated in the last islands of 
deconfined matter, as visualized in Fig. 

This baryon concentration may cause two interesting effects intensively debated 
during the last decade: either the quark islands stabilize themselves and survive 
as quark nugget relics or they convert into isothermal baryon density fluctuations. 
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Fig. 15 Left column: Growing bubbles of hadronic matter immersed in quark-gluon plasma. 
Right column: Shrinking islands of quark-plasma immersed in hadronic matter. For T <C 200 
MeV, the baryon charge is accumulated in the quark islands. 



We will discuss the former effect in the next subsection. The latter effect can af- 
fect the primordial nucleosynthesis if its survives up to then. A necessary condi- 
tion is / ^ 1 mlr^ioo A/ey , where / is the comoving scale characterizing the mean 
distance of the inhomogeneities, which is determined by the nucleation process via 
I ~ 0.8(Tc/100MeV) d„uc || with 

2 



13 



/145MeV\' 



4B 
T 



3/2 



cm: 



(41) 



V ) \ L ) V0.035B3/4. 

L is the latent heat in the bag model cf. left upper panel in Fig. |ll|. 

In confined matter the neutrons can rapidly diffuse thus causing an inhomogeneous 
neutron-to-proton ratio. Recent re-analyses Q let this scenario appear less prob- 
able: the low surface tension, discussed in subsection 6T, gives rise to a typical scale 
I ~ 1 cm and the detailed study of the nucleosynthesis shows no deviations from the 
homogeneous standard scenario. 

It should be mentioned, however, that the estimates quoted above rely on the 
homogeneous nucleation theory. Nucleation by impurities is more efficient and, as 
shown in p3| , can indeed render the inhomogeneity scale to / ~ 1 m, supposed there 
are pre-fabricated fluctuations or impurities on an appropriate scale, such as vortices, 
domain walls, magnetic monopoles, cosmic strings, or relic fluctuations from the elec- 
troweak transition. The required scale of Z ~ 1 m corresponds just to the red-shifted 
horizon scale at the electroweak transition (see below). In it is pointed out that 
temperature fluctuations, compatible with the COBE measurements, can give rise to 
a inhomogeneous nucleation with interesting consequences. 

Another idea is that a sufficiently large baryon concentration can survive up to 
now in form of quark nuggets being thus candidates for dark matter. For a recent 
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discussion cf. The quark nuggets must have baryon charges of ^ j^q40- -45 

because otherwise the baryon contrast would have diffused away. 

6.4 Strange matter 

Strangeness may play a particular role in the cosmic evolution. While nowadays the 
hadrons are mainly nucleons with the valence quark structure p — (uud), n — (ddu), 
just after confinement also pions and kaons were present. Kaons carry strangeness, 
e.g. = {us), but once the universe cools below fOO MeV they rapidly decay, and 
the strange particles disappear from the world. They can be afterwards newly created 
for a short while by strong interaction of high-energy cosmic rays with matter and in 
the laboratory. But the earlier universe was quite strange. Due to the long expansion 
and cooling time scales in comparison with the time scale of the weak interaction, 
there was chemical equilibrium according to the reactions changing strangeness 

d u + + De, s <-> u + + i^e- (42) 

This results in the balance equations 

= fJ-u + Me- - Mi^e ■ Ms = Mu + Me- " Mi'e : (43) 

from which fig — fid follows, or Us ~ rid- Thus, a substantial part of the quarks were 
strange ones (for a more detailed consideration including charge neutrality cf. p2| ). 
As pointed in ^ , within the bag model the strange quarks can stabilize the quark 
matter. This idea is re-analyzed with increasing sophistication and there seems to 
be a region in the parameter space still leaving the possibility of stable strange quark 
nuggets. These represent candidates for dark matter. 

The dedicated experiments with heavy- ion collisions, however, did not show posi- 
tive signals of strangelets . This constrains the parameter space further and makes 
the strangelet idea less probable. It should be stressed, however, that in heavy-ion 
collisions the time scales are much shorter and the conditions for creating strangelets 
may be less favorable. 

6.5 Effects of vanishing sound velocity 

As seen in Fig. |ll| the sound velocity vanishes in first-order phase transitions in the 
mixed phase. This observation has an important effect on the evolution of fluctuations: 
the restoring pressure gradient of density-enhanced regions vanishes, and these regions 
can collapse in free fall. 

6.5.1 Intermediate scales 

On intermediate scales, Z <C A ^ Rh, the fate of the fluctuations is followed by 
the cosmological perturbation formalism p[ . In particular, as pointed out in 
kinetically decoupled cold dark matter can be trapped in gravitational wells. Candi- 
dates for such dark matter are WIMPzillas^ {M > 10^ GeV), primordial black holes 
(10"^* < Af < Mq), and axions (M 10"^ eV) which can be accumulated 

^This nickname refers to the notion of Weakly Interacting Massive Particles, see |52| . 



B. Kampfer, Cosmic Phase Transitions 



23 



to lumps with mass of M ~ 1O^'^° "^°-'M0. Also primordial gravitational waves are 
modified, however, on scales today not accessible via pulsar timing. Afterwards, the 
neutrino diffusion is efficient enough to damp away the created density fluctuations 
till nucleosynthesis, unless the decoupled matter lumps. 

6.5.2 Horizon scales 

As pointed out in Q, on large scales, A ~ Rh, already previous causally non- 
connected weak fluctuations on super-horizon scales can, after entry into the horizon, 
collapse to black holes with M ^ IMq due to = during confinement (cf. for a 
further discussion). Such a mechanism can produce MACHO^ visible by gravitational 
lensing. 

6.5.3 Vanishing sound velocity during conGnement? 

Since both effects are related to the vanishing of the sound velocity, one has to analyze 
whether this may really happen. In the case of a strong crossover, instead of first- 
order phase transition, the sound velocity decreases but does not vanish completely. 
As shown in under such circumstances the resulting gravitational wells are too 
weak to substantially capture dark matter. 

One has also to mention that cj. — holds only for p — const in the phase mix- 
ture region when all chemical potentials fia vanish or only one independent chemical 
potential occurs, i.e., a = 1. In cosmic matter there are numerous potentials, such as 

Imposing (3 equilibrium, according to ( ^ ) and n p + e~De, and local electric 
charge neutrality, at least one of the chemical potentials changes discontinuously in the 
case of a naively constructed phase transition. Instead, as stressed in constructing 
the phase equilibrium with the averaged charge neutrality Uhad + ngcp = results 
(cf. Fig. in p 7^ const in the mixed state, and therefore the sound velocity does not 
strictly vanish. The effect of the various chemical potential may be small if /Xq, <C T 
and the sound velocity at least substantially drops in the mixed state. 

Otherwise one should notice that upon approaching the confinement temperature 
the sound velocity in the deconfined state itself already drops (cf. Fig. |l^). This is 
usually related to a softening of the equation of state. In ||5^ evidences of such a 
softening is extracted from data of heavy-ion collisions at AGS energies. It is matter 
of debate whether this is an indication of deconfinement effects at such comparatively 
low beam energies of 10 • • • 15 A-GeV. 

Moreover, the annihilation process of certain particle species is also related to a 
decrease of the sound velocity, cf. jsij. This turns out to be a marginal effect. To look 
into details let us consider the era when the muons disappear. At T ^ and at T <C 

the degeneracies are dcff — 14.25 and 10.75, respectively. By evaluating numerically 
the statistical integral (p5|), one gets the pressure as displayed in Fig. ^ Indeed, there 
is a small change of the scaled total pressure at T/m^ ^ 0.4. At T/m^ < 0.2, the 

muons can be neglected and at T/m^^ > 1 the radiation approximation, = 
with = 4, is sufficient. 

^ MACHO is an acronym for Massive Compact Halo Objects, see [p4]. 
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Fig. 16 Constructing the phase transition in presence of a some chemical potentials. Impos- 
ing charge neutrality yields the two curves labeled by QG and H for deconfined and confined 
matter, respectively. Pressure equality at fixed temperature and at equal neutron-baryon 
chemical potential /i„ (indicated by the dotted line) is connected with a jump of the electron 
chemical potential fie- Instead, averaged charge neutrality results in the gray S-shaped line. 
On this line, the pressures of both phases are equal but they change when moving along this 
line. 
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Fig. 17 The sound velocity of gluon matter as a function of temperature. Symbols depict 
the results of lattice QCD data from the group Bielefeld and the curve results from the 
quasi-particle model 
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Fig. 18 The scaled pressure as a function of the temperature of the muons (lower solid line: 
exact evaluation of the statistical integral ( p4| ) , dashed lines: using the first two terms of the 
high and low temperature expansion). Upper gray curve: total pressure taking into account 
the background of electrons, photons, and three neutrino species. 

Solving the evolution equations ( jl^ , p^ ) results in the temperature history as dis- 
played in the left panel of Fig. There is a very moderate change of the cooling due 
to muon annihilation. Only when incorrectly approximating the muon pressure by the 
first two terms in the high-temperature expansions one gets a more pronounced effect 
which looks like a phase transition (see right panel of Fig. pj| ). 

Thus the effect of particle annihilation is weak and unlikely to cause a reheating or 
a substantial softening of the equation of state. This can also be seen from the measure 
3p/e which drops from 1 to 0.95 during muon annihilation and rises to 1 afterwards. 
Therefore, particle annihilation can hardly cause strong effects, as mentioned in 

6.6 Interim resume 

The confinement transition is quite interesting, but the cosmic evolution is already slow 
compared with the strong interaction time scale, and probably the most severe effects 
are washed out. Quasi-equilibrium means memory loss. In addition some details are 
poorly known, such as the correct order of the transition and the surface tension (if 
any). Therefore, one might argue that earlier stages with more rapid evolution are out 
of equilibrium and do not leave such fragile imprints. 

7 Electroweak symmetry breaking 

In the electroweak standard model a la Glashow-Salam- Weinberg the Higgs field is 
responsible for the dynamical mass generation via spontaneous symmetry breaking. 
At sufficiently high temperatures, T > 100 GeV, the expectation value of the Higgs 
field is zero, i.e., the symmetry is restored and particles are massless. At T < 100 GeV 
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Fig. 19 The evolution of the temperature as a function of the scaled time. Lower dashed 
line: only the background, upper dashed line: background plus massless muons, solid line: 
background plus muon contribution. The left panel uses the correct evaluation of the pressure, 
while the right panel employs only the expansion described in text. 



the symmetry becomes broken because (Higgs) 7^ and particle masses become finite. 
For Higgs masses uih < 75 GeV, the symmetry breaking appears as weak first-order 
phase transition, while for the more hkely case of mn > 75 GeV a phase transition of 
second order happen s p7 , 58 1. Indeed, approaching mn 75 GeV, the surface tension 
drops dramatically ||5§[|. In view of the experience with the confinement transition, 
one may not expect drastic effects from this transition. 

However, there is an important issue, outside the scope of the present work, related 
to non-perturbative effects of sofitary field configurations called sphalerons: baryon 
and lepton numbers are not conserved. Therefore, at the electroweak transition the 
baryon and lepton excess of the universe might be generated. Since a final consistent 
picture does not seem to be elaborated yet, we refer to some recent surveys on this 
topic including also energies of the GUT scale, see [59t|. 



8 Inflation 

Equation (^) shows that at p < — the expansion of the universe becomes accel- 
erated. In particular, if the vacuum energy density dominates, i.e. e™'^ ^ Cq, and 
e = e™'^ = —p^^^ Ri const then R cx e''^*^™"*. A series of problems in modern cos- 
mology is solved by inflation if the scale factor R grows in the inflationary period by 
a factor of about 10^°. These standard problems, like horizon, flatness, homogeneity, 
age, and monopole problems, are summarized in many surveys, sec |ll6| . 

Inflation can be realized within a model with a dominating spatially homogeneous 
scalar field which, via a potential V{(f>), is self-interacting. The corresponding equation 
of state can be deduced from 

e = +F(0,r) + i02^ (44) 
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p = -F(0,T) + -02 (45) 

(quantum corrections cause a temperature dependence of the effective potential). For 
a sufficiently slow evolution, the kinetic term can be neglected and the coherent field 
mimics a perfect vacuum. Inflation is expected to happen on the GUT scale of T ~ 10^^ 
GeV and provides, via adiabatic Gaussian fluctuations, the seed for later structure 
formation. The various inflationary scenarios differ essentially only in the choice of 
the potential. To recall a few scenarios (cf. pOt): 

• old inflation: decay of a metastable ("wrong") vacuum caused by a potential 
with two local minima corresponding to a phase transition of first order, 

• new inflation: a potential which is very flat around (t>i = 0, where (pi is the initial 
state, and which has a minimum at (pji > 0, 

• chaotic inflation: a simple U-shaped potential and 0^ ^ evolving to (j^n = 0, 

• hybrid inflation: two coupled scalar fields with ip ^ (p, and 4> evolves in the slow 
rolling down potential of new inflation. 

We do not consider here scenarios beyond Einstein's theory or involving supersymme- 
try or string theory. 

8.1 Viscous inflation 

The original version of the old inflationary scenario ||6l|| assumed a supercooling by 30 
orders of magnitude to achieve the stretching of R by 30 order of magnitudes. Then, 
of course, the exponential grow of the space prevents a graceful exit to a homogeneous 
state after reheating. However, there are also phenomenological ansatze which realize 
successful inflation without such tremendous supercooling. For instance, in viscous 
inflation |62[ a suitably tuned volume viscosity can compensate the cooling, moreover, 
due to the replacement of the pure thermodynamic pressure p hy p — 3-^^ (see (^), 
the space expansion can accelerate for suitable values of ^. Then there happens a 
tremendous entropy production. To be specific, an increase of the comoving entropy, 
sR^, by a factor I0^° is sufficient to solve most of the standard problems. The key 
here is entropy production, which spoils the comoving entropy conservation and such 
an information link of very early and late stages of the universe. 

8.2 Tepid inflation 

Another way of entropy production is a phase transition with supercooling. Let us 
consider as an example an isothermal phase transition at T = const, as visualized in 
the left panel of Fig. |2^. The state of matter changes as indicated in the middle panel 
by the solid vertical arrow, i.e., there is a transition where the new phase (2) has larger 
pressure than the old phase (1). This causes a compression of the regions filled with 
the old phase. Taking into account the averaging prescription for the energy density 
during the phase transition 



e = cix + 62(1 — x) 



(46) 
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Fig. 20 Left panel: Isothermal phase transition with supercooling AT (heavy solid curve). 
Middle panel: Sketch of the equation of state near the equilibrium point "C". Right panel: 
Entropy increase as a function of supercooling. 



and analog equations for the pressure and the entropy density, the entropy increase 
for the perfect 2-phase fluid is found from (^ as 



{sR'j = i?'|^r|^(Pi-P2)i + (ri-T2) 



A + 3|a 



(47) 



A = XSi — (1 — x)s2, 

instead of (sR^)' — for the ideal 1-component fluid. To quantify the entropy produc- 
tion let us define an entropy increase factor by ^ = j-^- This factor is displayed in the 
right panel of Fig. ^ as a function of the supercooling. One observes that moderate 
supercooling in the order of 90 - 95% is sufficient for successful infiation. 

In | ]63| we have discussed a more dynamical model with a phase transition and 
an ansatz for the relaxation time approximation for X2{t). We find a sufficient tepid 
inflation. Now we consider the nucleation theory, presented in section |^, and study 
the evolution of the scale factor and the entropy increase. The results of a numerical 



integration of the evolution equations (|lj, |3^, |39|) are presented in Fig. |2lJ. There 
are two regimes: 

(i) a too small surface tension parameter ctq causes a too small supercooling and 
insufficient inflation, while 

(ii) a too large value of aa repeats the graceful exit problem of old inflation, i.e. the 
space is exponentially growing, and the growth of bubbles is too slow to fill the space, 
as visualized in Fig. p2[ 

However there is a small corridor in parameter space where we find power law 
inflation with R t"' , a > 1, resembling the extended inflation in Brans-Dicke theory. 
Numerically it is difficult to follow a successful inflationary evolution due to the needed 
fine tuning of ctq and the accuracy of the integration procedure with respect to the 
large number problems. Therefore, it is hard to make statements on the emerging 
fluctuation spectrum. 
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Fig. 21 Evolution of the scale factor (upper panel) and the entropy increase factor (lower 
panel) as a function of time. The dotted sections in the upper panel correspond to the flat 
sections in the lower panel and indicate that the old phase disappeared. Curves labeled by 
1, 2, 3, 4 correspond to ao = 2.064, 2.06499, 2.065, 2.066, respectively. For more details, like 
the parameterization of the equation of state, consult [p2, M. 




Fig. 22 Visualization of the graceful exit problem: the space expands too rapidly and the 
growing bubbles fail to percolate since nucleation is, due to large supercooling, too rare to 
fill the space. 
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9 Summary 

We should be aware that the evolution of the universe is accompanied by a permanent 
change of the states of matter, which can be in some interesting cases be considered as 
phase transitions. There is a change of certain order parameters, like a hypothetical 
scalar field (p on the GUT scale, the Higgs field on the electroweak scale, the chiral 
condensate on the QCD scale, or even now the quintessence. This causes partially 
drastic changes of matter like the transition from freely roaming quarks and gluons to 
hadrons, or the binding of nucleons into light nuclei during nucleosynthesis, or simply 
the annihilation of particles like e+e" — > 27. These changes induce also changes of the 
expansion dynamics of the universe, most drastic during inflation, e.g.. Furthermore 
in the transition from a radiation dominated universe to a matter dominated one, the 
time-evolution of the scale factor R{t) changes from 

uncertainty is introduced by the evidence recently found of a dominating vacuum 
energy density, the quintessence, because its back-extrapolation needs more definite 
knowledge on its nature. 

Assuming that in early stages quintessence is not dominating the energy density, 
we considered in some detail the confinement transition in the Big Bang and, very 
briefly, contrasted it to the Little Bang, presently investigated in heavy-ion collision 
experiments. In the Little Bang one observes direct electromagnetic radiation like real 
photons and virtual photons via 7* e+e". This points to temperature scales of 170 
MeV and larger, i.e., above the deconfinement transition. In contrast to this, there do 
not seem to be specific definite relics from the confinement era of the Big Bang, apart 
the protons and neutrons created in hadrosynthesis. We summarized possible relics, 
whose existence is hypothetical and the predictions of which are sometimes related to 
poorly known details of the confinement dynamics. 

Whereas the results of the nucleosynthesis processes specifically depend on the 
interplay of neutrino decoupling, neutron life time, and cosmic cooling rate, the 
hadrosynthesis at confinement seems to be rather unspecific. Nucleons emerge as 
other hadrons do as well, but the unstable hadrons decay shortly after confinement 
and the proton-to-neutron ratio is determined by /3 equilibrium. Hence, hadrosyn- 
thesis does not cause a specific hadron composition due to the cosmic evolution since 
chemical equilibrium means memory loss. 

Strictly speaking, the success of the primordial nucleosynthesis makes us strongly 
believe in a maximum temperature of Tn > 1 MeV, while the prediction of a flat 
universe and of an appropriate fluctuation spectrum for structure formation by infla- 
tionary scenarios makes us believe in Tn > 10"'^^ GeV. Within these scales, a consistent 
picture of the evolution of the universe seems to emerge, but it leaves a consistent treat- 
ment of the "initial state" and the reheating (or, better, "pre-heating" ) to Tr as the 
most challenging problem probably to be solved within quantum gravity. 
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